The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome.

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

Integration of beef cattle international pedigree and genomic estimated breeding values into national evaluations, with an application to the Italian Limousin population.

BACKGROUND: International evaluations combine data from different countries allowing breeders to have access to larger panels of elite bulls and to increase the accuracy of estimated breeding values (EBV). However, international and national evaluations can use different sources of information to compute EBV (EBVINT and EBVNAT, respectively), leading to differences between them. Choosing one of these EBV results in losing the information that is contained only in the discarded EBV. Our objectives were to define and validate a procedure to integrate publishable sires' EBVINT and their associated reliabilities computed from pedigree-based or single-step international beef cattle evaluations into national evaluations to obtain "blended" EBV. The Italian (ITA) pedigree-based national evaluation was used as a case study to validate the integration procedure. METHODS: Publishable sires' international information, i.e. EBVINT and their associated reliabilities, was included in the national evaluation as pseudo-records. Data were available for 444,199 individual age-adjusted weaning weights of Limousin cattle from eight countries and 17,607 genotypes from four countries (ITA excluded). To mimic differences between international and national evaluations, international evaluations included phenotypes (and genotypes) of animals born prior to January 2019, while national evaluations included ITA phenotypes of animals born until April 2019. International evaluations using all available information were considered as reference scenarios. Publishable sires were divided into three groups: sires with ≥ 15, < 15 and no recorded offspring in ITA. RESULTS: Overall, for these three groups, integrating either pedigree-based or single-step international information into national pedigree-based evaluations improved the similarity of the blended EBV with the reference EBV compared to national evaluations without integration. For instance, the correlation with the reference EBV for direct (maternal) EBV went from 0.61 (0.79) for a national evaluation without integration to 0.97 (0.88) when integrating single-step international information, on average across all groups of publishable sires. CONCLUSIONS: Our proposed one-animal-at-a-time integration procedure yields blended EBV that are in close agreement with full international EBV for all groups of animals analysed. The procedure can be directly applied by countries since it does not rely on specific software and is computationally inexpensive, allowing straightforward integration of publishable sires' EBVINT from pedigree-based or single-step based international beef cattle evaluations into national evaluations.

High-resolution structural variants catalogue in a large-scale whole genome sequenced bovine family cohort data.

BACKGROUND: Structural variants (SVs) are chromosomal segments that differ between genomes, such as deletions, duplications, insertions, inversions and translocations. The genomics revolution enabled the discovery of sub-microscopic SVs via array and whole-genome sequencing (WGS) data, paving the way to unravel the functional impact of SVs. Recent human expression QTL mapping studies demonstrated that SVs play a disproportionally large role in altering gene expression, underlining the importance of including SVs in genetic analyses. Therefore, this study aimed to generate and explore a high-quality bovine SV catalogue exploiting a unique cattle family cohort data (total 266 samples, forming 127 trios). RESULTS: We curated 13,731 SVs segregating in the population, consisting of 12,201 deletions, 1,509 duplications, and 21 multi-allelic CNVs (> 50-bp). Of these, we validated a subset of copy number variants (CNVs) utilising a direct genotyping approach in an independent cohort, indicating that at least 62% of the CNVs are true variants, segregating in the population. Among gene-disrupting SVs, we prioritised two likely high impact duplications, encompassing ORM1 and POPDC3 genes, respectively. Liver expression QTL mapping results revealed that these duplications are likely causing altered gene expression, confirming the functional importance of SVs. Although most of the accurately genotyped CNVs are tagged by single nucleotide polymorphisms (SNPs) ascertained in WGS data, most CNVs were not captured by individual SNPs obtained from a 50K genotyping array. CONCLUSION: We generated a high-quality SV catalogue exploiting unique whole genome sequenced bovine family cohort data. Two high impact duplications upregulating the ORM1 and POPDC3 are putative candidates for postpartum feed intake and hoof health traits, thus warranting further investigation. Generally, CNVs were in low LD with SNPs on the 50K array. Hence, it remains crucial to incorporate CNVs via means other than tagging SNPs, such as investigation of tagging haplotypes, direct imputation of CNVs, or direct genotyping as done in the current study. The SV catalogue and the custom genotyping array generated in the current study will serve as valuable resources accelerating utilisation of full spectrum of genetic variants in bovine genomes.

Classifying aneuploidy in genotype intensity data using deep learning.

Aneuploidy is the loss or gain of one or more chromosomes. Although it is a rare phenomenon in liveborn individuals, it is observed in livestock breeding populations. These breeding populations are often routinely genotyped and the genotype intensity data from single nucleotide polymorphism (SNP) arrays can be exploited to identify aneuploidy cases. This identification is a time-consuming and costly task, because it is often performed by visual inspection of the data per chromosome, usually done in plots of the intensity data by an expert. Therefore, we wanted to explore the feasibility of automated image classification to replace (part of) the visual detection procedure for any diploid species. The aim of this study was to develop a deep learning Convolutional Neural Network (CNN) classification model based on chromosome level plots of SNP array intensity data that can classify the images into disomic, monosomic and trisomic cases. A multispecies dataset enriched for aneuploidy cases was collected containing genotype intensity data of 3321 disomic, 1759 monosomic and 164 trisomic chromosomes. The final CNN model had an accuracy of 99.9%, overall precision was 1, recall was 0.98 and the F1 score was 0.99 for classifying images from intensity data. The high precision assures that cases detected are most likely true cases, however, some trisomy cases may be missed (the recall of the class trisomic was 0.94). This supervised CNN model performed much better than an unsupervised k-means clustering, which reached an accuracy of 0.73 and had especially difficult to classify trisomic cases correctly. The developed CNN classification model provides high accuracy to classify aneuploidy cases based on images of plotted X and Y genotype intensity values. The classification model can be used as a tool for routine screening in large diploid populations that are genotyped to get a better understanding of the incidence and inheritance, and in addition, avoid anomalies in breeding candidates.

Impact of genomic preselection on subsequent ssGBLUP evaluation of preselected animals for scarcely recorded feed intake in pigs.

We have previously shown that single-step genomic best linear unbiased prediction (ssGBLUP) estimates breeding values of genomically preselected animals without preselection bias for widely recorded traits, that is traits recorded for the majority of animals in the breeding population. This study investigated the impact of genomic preselection (GPS) on accuracy and bias in ssGBLUP evaluation of genomically preselected animals for a scarcely recorded trait, that is a trait recorded for only a small proportion of the animals, which generally has a lower prediction accuracy than widely recorded traits, mainly due to having a much smaller number of phenotypes available. We used data from a commercial pig breeding program, considering feed intake as a scarcely recorded target trait, being available for ~30% of the animals with phenotypes for any trait, and average daily gain, backfat thickness and loin depth as widely recorded predictor traits, being available for >95% of the animals with phenotypes for any trait. The data contained the routine GPS implemented by commercial animal breeding programs, and we retrospectively implemented two scenarios with additional layers of GPS by discarding pedigree, genotypes and phenotypes of animals without progeny. The ssGBLUP evaluation following GPS used records only from the target trait, only from the predictor traits, or both. Accuracy for feed intake did not differ statistically across GPS scenarios, although it tended to decrease with more intense GPS. The accuracy had average values of 0.37, 0.44, and 0.45 across all GPS scenarios when, respectively, records from only the target trait, only the predictor traits, or both were used in the ssGBLUP evaluation. Considerable deflation of the genomic breeding values for feed intake was observed in the most stringent GPS scenario, due to the variance components being underestimated as a result of the limited amount of strongly preselected data. As long as (co)variance components were unbiased, no or only marginal bias was observed. These results for accuracy and bias were observed whether records of the scarcely recorded target trait, of the predictor traits, or both were used in the ssGBLUP evaluation. Our results show that for the scarcely recorded feed intake in pigs, ssGBLUP is able to estimate breeding values of preselected animals without preselection bias, similarly as previously observed for widely recorded traits.

International single-step SNPBLUP beef cattle evaluations for Limousin weaning weight.

BACKGROUND: Compared to national evaluations, international collaboration projects further improve accuracies of estimated breeding values (EBV) by building larger reference populations or performing a joint evaluation using data (or proxy of them) from different countries. Genomic selection is increasingly adopted in beef cattle, but, to date, the benefits of including genomic information in international evaluations have not been explored. Our objective was to develop an international beef cattle single-step genomic evaluation and investigate its impact on the accuracy and bias of genomic evaluations compared to current pedigree-based evaluations. METHODS: Weaning weight records were available for 331,593 animals from seven European countries. The pedigree included 519,740 animals. After imputation and quality control, 17,607 genotypes at a density of 57,899 single nucleotide polymorphisms (SNPs) from four countries were available. We implemented two international scenarios where countries were modelled as different correlated traits: an international genomic single-step SNP best linear unbiased prediction (SNPBLUP) evaluation (ssSNPBLUPINT) and an international pedigree-based BLUP evaluation (PBLUPINT). Two national scenarios were implemented for pedigree and genomic evaluations using only nationally submitted phenotypes and genotypes. Accuracies, level and dispersion bias of EBV of animals born from 2014 onwards, and increases in population accuracies were estimated using the linear regression method. RESULTS: On average across countries, 39 and 17% of sires and maternal-grand-sires with recorded (grand-)offspring across two countries were genotyped. ssSNPBLUPINT showed the highest accuracies of EBV and, compared to PBLUPINT, led to increases in population accuracy of 13.7% for direct EBV, and 25.8% for maternal EBV, on average across countries. Increases in population accuracies when moving from national scenarios to ssSNPBLUPINT were observed for all countries. Overall, ssSNPBLUPINT level and dispersion bias remained similar or slightly reduced compared to PBLUPINT and national scenarios. CONCLUSIONS: International single-step SNPBLUP evaluations are feasible and lead to higher population accuracies for both large and small countries compared to current international pedigree-based evaluations and national evaluations. These results are likely related to the larger multi-country reference population and the inclusion of phenotypes from relatives recorded in other countries via single-step international evaluations. The proposed international single-step approach can be applied to other traits and breeds.

Impact of genomic preselection on subsequent genetic evaluations with ssGBLUP using real data from pigs.

BACKGROUND: Empirically assessing the impact of preselection on genetic evaluation of preselected animals requires comparing scenarios that take different approaches into account, including scenarios without preselection. However, preselection is almost always performed in animal breeding programs, so it is difficult to have a dataset without preselection. Hence, most studies on preselection have used simulated datasets, and have concluded that genomic estimated breeding values (GEBV) from subsequent single-step genomic best linear unbiased prediction (ssGBLUP) evaluations are unbiased. The aim of this study was to investigate the impact of genomic preselection (GPS) on accuracy and bias in subsequent ssGBLUP evaluations, using data from a commercial pig breeding program. METHODS: We used data on average daily gain during performance testing, average daily gain throughout life, backfat thickness, and loin depth from one sire line and one dam line of pigs. As these traits have different weights in the breeding goals of the two lines, we analyzed the lines separately. For each line, we implemented a reference GPS scenario that kept all available data, against which the next two scenarios were compared. We then implemented two other scenarios with additional layers of GPS by removing all animals without progeny either (i) only in the validation generation, or (ii) in all generations. We conducted subsequent ssGBLUP evaluations for each GPS scenario, using all the data remaining after implementing the GPS scenario. Accuracy and bias were computed by comparing GEBV against progeny yield deviations of validation animals. RESULTS: Results for all traits and in both lines showed a marginal loss in accuracy due to the additional layers of GPS. Average accuracies across all GPS scenarios in the two lines were 0.39, 0.47, 0.56, and 0.60, for average daily gain during performance testing and throughout life, backfat thickness, and loin depth, respectively. Biases were largely absent, and when present, did not differ greatly between the GPS scenarios. CONCLUSIONS: We conclude that the impact of preselection on accuracy and bias in subsequent ssGBLUP evaluations of selection candidates in pigs is generally minimal. We expect this conclusion to apply for other animal breeding programs as well, since preselection of any type or intensity generally has the same effect in animal breeding programs.

Development of resilience indicator traits based on daily step count data for dairy cattle breeding.

BACKGROUND: Resilient animals are minimally affected by disturbances, such as diseases and heat stress, and quickly recover. Daily activity data can potentially indicate resilience, because resilient animals likely keep variations due to disturbances that threat animal homeostasis at a low magnitude. We used daily step count of cows to define resilience indicators based on theory, exploratory analysis and literature, and then investigated if they can be used to genetically improve resilience by estimating heritability and repeatability, and genetic associations with other resilience-related traits, i.e. health traits, longevity, fertility, and body condition score (BCS). RESULTS: Two groups of resilience indicators were defined: indicators describing (1) mean step count at different lactation stages for individual cows, and (2) fluctuations in step count from individual step count curves. Heritability estimates were highest for resilience indicators describing mean step count, from 0.22 for the 2-week period pre-partum to 0.45 for the whole lactation. High mean step count was consistently, but weakly, genetically correlated with good health, fertility, and longevity, and high BCS. Heritability estimates of resilience indicators describing fluctuations ranged from 0.01 for number of step count drops to 0.15 for the mean of negative residuals from individual curves. Genetic correlations with health traits, longevity, fertility, and BCS were mostly weak, but were moderate and favorable for autocorrelation of residuals (- 0.33 to - 0.44) and number of step count drops (- 0.44 to - 0.56) with hoof health, fertility, and BCS. Resilience indicators describing variability of residuals and mean of negative residuals showed strong genetic correlations with mean step count (0.86 to 0.95, absolute), which suggests that adjustment for step count level is needed. After adjustment, 'mean of negative residuals' was highly genetically correlated with hoof health, fertility, and BCS. CONCLUSIONS: Mean step count, autocorrelation and mean of negative residuals showed most potential as resilience indicators based on resilience theory, heritability, and genetic associations with health, fertility, and body condition score. Other resilience indicators were heritable, but had unfavorable genetic correlations with several health traits. This study is an important first step in the exploration of the use of activity data to breed more resilient livestock.

The impact of direct-maternal genetic correlations on international beef cattle evaluations for Limousin weaning weight.

In beef cattle maternally influenced traits, estimates of direct-maternal genetic correlations (rdm) are usually reported to be negative. In international evaluations, rdm can differ both within countries (rdm_WC) and between countries (rdm_BC). The rdm_BC are difficult to estimate and are assumed to be zero in the current model for international beef cattle evaluations (Interbeef). Our objective was to investigate re-ranking of international estimated breeding values (IEBVs) in international beef cattle evaluations between models that either used estimated values for rdm or assumed them to be 0. Age-adjusted weaning weights and pedigree data were available for Limousin beef cattle from ten European countries. International EBVs were obtained using a multi-trait animal model with countries modeled as different traits. We compared IEBVs from a model that uses estimated rdm_BC (ranging between -0.14 and +0.14) and rdm_WC (between -0.33 and +0.40) with IEBVs obtained either from the current model that assumes rdm_BC to be 0, or from an alternative model that assumes both rdm_BC and rdm_WC to be 0. Direct and maternal IEBVs were compared across those three scenarios for different groups of animals. The ratio of population accuracies from the linear regression method was used to further investigate the impact of rdm on international evaluations, for both the whole set of animals in the evaluation and the domestic ones. Ignoring rdm_BC, i.e., replacing estimated values with 0, resulted in no (rank correlations > 0.99) or limited (between 0.98 and 0.99) re-ranking for direct and maternal IEBVs, respectively. Both rdm_BC and rdm_WC had less impact on direct IEBVs than on maternal IEBVs. Re-ranking of maternal IEBVs decreased with increasing reliability. Ignoring rdm_BC resulted in no re-ranking for sires with IEBVs that might be exchanged across countries and limited re-ranking for the top 100 sires. Using estimated rdm_BC values instead of considering them to be 0 resulted in null to limited increases in population accuracy. Ignoring both rdm_BC and rdm_WC resulted in considerable re-ranking of animals' IEBVs in all groups of animals evaluated. This study showed the limited impact of the current practice of ignoring rdm_BC in international evaluations for Limousin weaning weight, most likely because the estimated rdm_BC was close to 0. We expect that these conclusions can be extended to other traits that have reported rdm values in the range of rdm_WC values for weaning weight in Limousin.

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle.

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.

Breeding goal traits accounting for feed intake capacity and roughage or concentrate intake separately.

Current breeding tools aiming to improve feed efficiency use definitions based on total dry matter intake (DMI); for example, residual feed intake or feed saved. This research aimed to define alternative traits using existing data that differentiate between feed intake capacity and roughage or concentrate intake, and to investigate the phenotypic and genetic relationships among these traits. The data set contained 39,017 weekly milk yield, live weight, and DMI records of 3,164 cows. The 4 defined traits were as follows: (1) Feed intake capacity (FIC), defined as the difference between how much a cow ate and how much she was expected to eat based on diet satiety value and status of the cow (parity and lactation stage); (2) feed saved (FS), defined as the difference between the measured and the predicted DMI, based on the regression of DMI on milk components within experiment; (3) residual roughage intake (RRI), defined as the difference between the measured and the predicted roughage intake, based on the regression of roughage intake on milk components and concentrate intake within experiment; and (4) residual concentrate intake (RCI), defined as the difference between the measured and the predicted concentrate intake, based on the regression of concentrate intake on milk components and roughage intake within experiment. The phenotypic correlations were -0.72 between FIC and FS, -0.84 between FS and RRI, and -0.53 between FS and RCI. Heritability of FIC, FS, RRI, and RCI were estimated to be 0.21, 0.12, 0.15, and 0.03, respectively. The genetic correlations were -0.81 between FS and FIC, -0.96 between FS and RRI, and -0.25 between FS and RCI. Concentrate intake and RCI had low heritability. Genetic correlation between DMI and FIC was 0.98. Although the defined traits had moderate phenotypic correlations, the genetic correlations between DMI, FS, FIC, and RRI were above 0.79 (in absolute terms), suggesting that these traits are genetically similar. Therefore, selecting for FIC is expected to simply increase DMI and RRI, and there seems to be little advantage in separating concentrate and roughage intake in the genetic evaluation, because measured concentrate intake was determined by the feeding system in our data and not by the genetics of the cow.

Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle.

Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10-5) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10-8) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.

Avoiding preselection bias in subsequent single-step genomic BLUP evaluations of genomically preselected animals.

In animal breeding, parents of the next generation are usually selected in multiple stages, and the initial stages of this selection are called preselection. Preselection reduces the information available for subsequent evaluation of preselected animals and this sometimes leads to bias. The objective of this study was to establish the minimum information required to subsequently evaluate genomically preselected animals without bias arising from preselection, with single-step genomic best linear unbiased prediction (ssGBLUP). We simulated a nucleus of a breeding program in which a recent population of 15 generations was produced. In each generation, parents of the next generation were selected in a single-stage selection based on pedigree BLUP. However, in generation 15, 10% of male and 15% of female offspring were preselected on their genomic estimated breeding values (GEBV). These GEBV were estimated using ssGBLUP, including the pedigree of all animals in generations 0-15, genotypes of all animals in generations 13-15 and phenotypes of all animals in generations 11-14. In subsequent ssGBLUP evaluation of these preselected animals, genotypes and phenotypes from various groups of animals were excluded one after another. We found that GEBV of the preselected animals were only estimated without preselection bias when genotypes and phenotypes of all animals in generations 13 and 14 and of the preselected animals were included in the subsequent evaluation. We also found that genotypes of the animals discarded at preselection only helped in reducing preselection bias in GEBV of their preselected sibs when genotypes of their parents were absent or excluded from the subsequent evaluation. We concluded that to prevent preselection bias in subsequent ssGBLUP evaluation of genomically preselected animals, information representative of the reference data used in the evaluation at preselection and genotypes and phenotypes of the preselected animals are needed in the subsequent evaluation.

Inbreeding depression across the genome of Dutch Holstein Friesian dairy cattle.

BACKGROUND: Inbreeding depression refers to the decrease in mean performance due to inbreeding. Inbreeding depression is caused by an increase in homozygosity and reduced expression of (on average) favourable dominance effects. Dominance effects and allele frequencies differ across loci, and consequently inbreeding depression is expected to differ along the genome. In this study, we investigated differences in inbreeding depression across the genome of Dutch Holstein Friesian cattle, by estimating dominance effects and effects of regions of homozygosity (ROH). METHODS: Genotype (75 k) and phenotype data of 38,792 cows were used. For nine yield, fertility and udder health traits, GREML models were run to estimate genome-wide inbreeding depression and estimate additive, dominance and ROH variance components. For this purpose, we introduced a ROH-based relationship matrix. Additive, dominance and ROH effects per SNP were obtained through back-solving. In addition, a single SNP GWAS was performed to identify significant additive, dominance or ROH associations. RESULTS: Genome-wide inbreeding depression was observed for all yield, fertility and udder health traits. For example, a 1% increase in genome-wide homozygosity was associated with a decrease in 305-d milk yield of approximately 99 kg. For yield traits only, including dominance and ROH effects in the GREML model resulted in a better fit (P < 0.05) than a model with only additive effects. After correcting for the effect of genome-wide homozygosity, dominance and ROH variance explained less than 1% of the phenotypic variance for all traits. Furthermore, dominance and ROH effects were distributed evenly along the genome. The most notable region with a favourable dominance effect for yield traits was on chromosome 5, but overall few regions with large favourable dominance effects and significant dominance associations were detected. No significant ROH-associations were found. CONCLUSIONS: Inbreeding depression was distributed quite equally along the genome and was well captured by genome-wide homozygosity. These findings suggest that, based on 75 k SNP data, there is little benefit of accounting for region-specific inbreeding depression in selection schemes.

Heritability estimates of the novel trait 'suppressed in ovo virus infection' in honey bees (Apis mellifera).

Honey bees are under pressure due to abnormal high colony death rates, especially during the winter. The infestation by the Varroa destructor mite and the viruses that this ectoparasite transmits are generally considered as the bees' most important biological threats. Almost all efforts to remedy this dual infection have so far focused on the control of the Varroa mite alone and not on the viruses it transmits. In the present study, the sanitary control of breeding queens was conducted on eggs taken from drone brood for 4 consecutive years (2015-2018). The screening was performed on the sideline of an ongoing breeding program, which allowed us to estimate the heritabilities of the virus status of the eggs. We used the term 'suppressed in ovo virus infection' (SOV) for this novel trait and found moderate heritabilities for the presence of several viruses simultaneously and for the presence of single viral species. Colonies that expressed the SOV trait seemed to be more resilient to virus infections as a whole with fewer and less severe Deformed wing virus infections in most developmental stages, especially in the male caste. The implementation of this novel trait into breeding programs is recommended.

Using prior information from humans to prioritize genes and gene-associated variants for complex traits in livestock.

Genome-Wide Association Studies (GWAS) in large human cohorts have identified thousands of loci associated with complex traits and diseases. For identifying the genes and gene-associated variants that underlie complex traits in livestock, especially where sample sizes are limiting, it may help to integrate the results of GWAS for equivalent traits in humans as prior information. In this study, we sought to investigate the usefulness of results from a GWAS on human height as prior information for identifying the genes and gene-associated variants that affect stature in cattle, using GWAS summary data on samples sizes of 700,000 and 58,265 for humans and cattle, respectively. Using Fisher's exact test, we observed a significant proportion of cattle stature-associated genes (30/77) that are also associated with human height (odds ratio = 5.1, p = 3.1e-10). Result of randomized sampling tests showed that cattle orthologs of human height-associated genes, hereafter referred to as candidate genes (C-genes), were more enriched for cattle stature GWAS signals than random samples of genes in the cattle genome (p = 0.01). Randomly sampled SNPs within the C-genes also tend to explain more genetic variance for cattle stature (up to 13.2%) than randomly sampled SNPs within random cattle genes (p = 0.09). The most significant SNPs from a cattle GWAS for stature within the C-genes did not explain more genetic variance for cattle stature than the most significant SNPs within random cattle genes (p = 0.87). Altogether, our findings support previous studies that suggest a similarity in the genetic regulation of height across mammalian species. However, with the availability of a powerful GWAS for stature that combined data from 8 cattle breeds, prior information from human-height GWAS does not seem to provide any additional benefit with respect to the identification of genes and gene-associated variants that affect stature in cattle.

Using short read sequencing to characterise balanced reciprocal translocations in pigs.

BACKGROUND: A balanced constitutional reciprocal translocation (RT) is a mutual exchange of terminal segments of two non-homologous chromosomes without any loss or gain of DNA in germline cells. Carriers of balanced RTs are viable individuals with no apparent phenotypical consequences. These animals produce, however, unbalanced gametes and show therefore reduced fertility and offspring with congenital abnormalities. This cytogenetic abnormality is usually detected using chromosome staining techniques. The aim of this study was to test the possibilities of using paired end short read sequencing for detection of balanced RTs in boars and investigate their breakpoints and junctions. RESULTS: Balanced RTs were recovered in a blinded analysis, using structural variant calling software DELLY, in 6 of the 7 carriers with 30 fold short read paired end sequencing. In 15 non-carriers we did not detect any RTs. Reducing the coverage to 20 fold, 15 fold and 10 fold showed that at least 20 fold coverage is required to obtain good results. One RT was not detected using the blind screening, however, a highly likely RT was discovered after unblinding. This RT was located in a repetitive region, showing the limitations of short read sequence data. The detailed analysis of the breakpoints and junctions suggested three junctions showing microhomology, three junctions with blunt-end ligation, and three micro-insertions at the breakpoint junctions. The RTs detected also showed to disrupt genes. CONCLUSIONS: We conclude that paired end short read sequence data can be used to detect and characterize balanced reciprocal translocations, if sequencing depth is at least 20 fold coverage. However, translocations in repetitive areas may require large fragments or even long read sequence data.

Investigating the impact of preselection on subsequent single-step genomic BLUP evaluation of preselected animals.

BACKGROUND: Preselection of candidates, hereafter referred to as preselection, is a common practice in breeding programs. Preselection can cause bias and accuracy loss in subsequent pedigree-based best linear unbiased prediction (PBLUP). However, the impact of preselection on subsequent single-step genomic BLUP (ssGBLUP) is not completely clear yet. Therefore, in this study, we investigated, across different heritabilities, the impact of intensity and type of preselection on subsequent ssGBLUP evaluation of preselected animals. METHODS: We simulated a nucleus of a breeding programme, in which a recent population of 15 generations was produced with PBLUP-based selection. In generation 15 of this recent population, the parents of the next generation were preselected using several preselection scenarios. These scenarios were combinations of three intensities of preselection (no, high or very high preselection) and three types of preselection (genomic, parental average or random), across three heritabilities (0.5, 0.3 or 0.1). Following each preselection scenario, a subsequent evaluation was performed using ssGBLUP by excluding all the information from the preculled animals, and these genetic evaluations were compared in terms of accuracy and bias for the preselected animals, and in terms of realized genetic gain. RESULTS: Type of preselection affected selection accuracy at both preselection and subsequent evaluation stages. While preselection accuracy decreased, accuracy in the subsequent ssGBLUP evaluation increased, from genomic to parent average to random preselection scenarios. Bias was always negligible. Genetic gain decreased from genomic to parent average to random preselection scenarios. Genetic gain also decreased with increasing intensity of preselection, but only by a maximum of 0.1 additive genetic standard deviation from no to very high genomic preselection scenarios. CONCLUSIONS: Using ssGBLUP in subsequent evaluations prevents preselection bias, irrespective of intensity and type of preselection, and heritability. With GPS, in addition to reducing the phenotyping effort considerably, the use of ssGBLUP in subsequent evaluations realizes only a slightly lower genetic gain than that realized without preselection. This is especially the case for traits that are expensive to measure (e.g. feed intake of individual broiler chickens), and traits for which phenotypes can only be measured at advanced stages of life (e.g. litter size in pigs).

Impact of sub-setting the data of the main Limousin beef cattle population on the estimates of across-country genetic correlations.

BACKGROUND: Cattle international genetic evaluations allow the comparison of estimated breeding values (EBV) across different environments, i.e. countries. For international evaluations, across-country genetic correlations (rg) need to be estimated. However, lack of convergence of the estimated parameters and high standard errors of the rg are often experienced for beef cattle populations due to limited across-country genetic connections. Furthermore, using all available genetic connections to estimate rg is prohibitive due to computational constraints, thus sub-setting the data is necessary. Our objective was to investigate and compare the impact of strategies of data sub-setting on estimated across-country rg and their computational requirements. METHODS: Phenotype and pedigree information for age-adjusted weaning weight was available for ten European countries and 3,128,338 Limousin beef cattle males and females. Using a Monte Carlo based expectation-maximization restricted maximum likelihood (MC EM REML) methodology, we estimated across-country rg by using a multi-trait animal model where countries are modelled as different correlated traits. Values of rg were estimated using the full data and four different sub-setting strategies that aimed at selecting the most connected herds from the largest population. RESULTS: Using all available data, direct and maternal rg (standard errors in parentheses) were on average equal to 0.79 (0.14) and 0.71 (0.19), respectively. Direct-maternal within-country and between-country rg were on average equal to - 0.12 (0.09) and 0.00 (0.14), respectively. Data sub-setting scenarios gave similar results: on average, estimated rg were smaller compared to using all data for direct (0.02) and maternal (0.05) genetic effects. The largest differences were obtained for the direct-maternal within-country and between-country rg, which were, on average 0.13 and 0.12 smaller compared to values obtained by using all data. Standard errors always increased when reducing the data, by 0.02 to 0.06, on average. The proposed sub-setting strategies reduced the required computing time up to 22% compared to using all data. CONCLUSIONS: Estimating all 120 across-country rg that are required for beef cattle international evaluations, using a multi-trait MC EM REML approach, is feasible but involves long computing time. We propose four strategies to reduce computational requirements while keeping a multi-trait estimation approach. In all scenarios with data sub-setting, the estimated rg were consistently smaller (mainly for direct-maternal rg) and had larger standard errors.

A deterministic equation to predict the accuracy of multi-population genomic prediction with multiple genomic relationship matrices.

BACKGROUND: A multi-population genomic prediction (GP) model in which important pre-selected single nucleotide polymorphisms (SNPs) are differentially weighted (MPMG) has been shown to result in better prediction accuracy than a multi-population, single genomic relationship matrix ([Formula: see text]) GP model (MPSG) in which all SNPs are weighted equally. Our objective was to underpin theoretically the advantages and limits of the MPMG model over the MPSG model, by deriving and validating a deterministic prediction equation for its accuracy. METHODS: Using selection index theory, we derived an equation to predict the accuracy of estimated total genomic values of selection candidates from population [Formula: see text] ([Formula: see text]), when individuals from two populations, [Formula: see text] and [Formula: see text], are combined in the training population and two [Formula: see text], made respectively from pre-selected and remaining SNPs, are fitted simultaneously in MPMG. We used simulations to validate the prediction equation in scenarios that differed in the level of genetic correlation between populations, heritability, and proportion of genetic variance explained by the pre-selected SNPs. Empirical accuracy of the MPMG model in each scenario was calculated and compared to the predicted accuracy from the equation. RESULTS: In general, the derived prediction equation resulted in accurate predictions of [Formula: see text] for the scenarios evaluated. Using the prediction equation, we showed that an important advantage of the MPMG model over the MPSG model is its ability to benefit from the small number of independent chromosome segments ([Formula: see text]) due to the pre-selected SNPs, both within and across populations, whereas for the MPSG model, there is only a single value for [Formula: see text], calculated based on all SNPs, which is very large. However, this advantage is dependent on the pre-selected SNPs that explain some proportion of the total genetic variance for the trait. CONCLUSIONS: We developed an equation that gives insight into why, and under which conditions the MPMG outperforms the MPSG model for GP. The equation can be used as a deterministic tool to assess the potential benefit of combining information from different populations, e.g., different breeds or lines for GP in livestock or plants, or different groups of people based on their ethnic background for prediction of disease risk scores.